Anderson Localization of cold atomic gases with effective spin-orbit interaction in a 

quasiperiodic optical lattice 
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We theoretically investigate the localization properties of a spin-orbit coupled spin- 1/2 particle 
moving in a one-dimensional quasiperiodic potential, which can be experimentally implemented us- 
ing cold atoms trapped in a quasiperiodic optical lattice potential and external laser fields. We 
present the phase diagram in the parameter space of the disorder strength and those related to the 
spin-orbit coupling. The phase diagram is verified via multifractal analysis of the atomic wavefunc- 
tions and the numerical simulation of diffusion dynamics. We found that spin-orbit coupling can 
lead to the spectra mixing (coexistence of extended and localized states) and the appearance of 
mobility edges. 
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I. INTRODUCTION 

Anderson localization (AL) is considered as a funda- 
mental physical phenomenon, which was first studied in 
the system of noninteracting electrons in a crystal with 
impurities fl*!. AL predicts the absence of diffusion of 
electronic spin, which stems from the disorder of crystal 
and is the result of quantum interference. Since disorder 
is ubiquitous, AL is rather universal and can occur in a 
variety of other physical systems including light waves [2] 
and atomic matter- waves Due to its intimate re- 

lation with metal-insulator transition, many interesting 
topics in AL such as the interplay between nonlinearity 
and disorder are still under intense study. 

In condensed-matter physics, spin-orbit (SO) coupling 
originates from the interaction between the intrinsic spin 
of an electron and the magnetic field induced by its move- 
ment. It connects the electronic spin to its orbital motion 
and thus the electron transport becomes spin-dependent. 
SO interaction can significantly affect AL and this prob- 
lem had been addressed by a few works in the electronic 
gas system [l^, [llj . 

The experimental realization of ultracold quantum 
gases, together with the technique of optical lattice po- 
tential, have provided a powerful playground for the sim- 
ulation of condensed-matter systems. In this compos- 
ite system, one can achieve unprecedented control over 
almost all parameters by optical or magnetical means. 
Specifically, pseudo-disorder can be generated by super- 
imposing two standing optical waves of incommensurate 
wavelengths together. As a consequence, AL of atomic 
matter wave can take place, which had been experimen- 
tally observed [1, S 111 and extensively studied in theory 

id- 

This work is motivated by the recent experimental re- 
alization of SO coupling in ultracold atomic gas 13- 16]. 
We investigate the impact of SO coupling on Anderson 
localization of a spin-1/2 particle using the system of cold 



atomic gases trapped in a quasiperiodic one-dimensional 
(ID) optical lattice potential and simultaneously subject 
to the laser-induced SO interaction. The similar topic 
had also been addressed in [13, [111 , with the focus on the 
localization properties of relativistic Dirac particles with 
cold atoms in a light-induced gauge field. Our model and 
method are different from theirs and we do not consider 
the relativistic region. 

The paper is organized as follows. Sec. [Il] introduces 
the theoretical model and tight-binding approximation 
is applied in Sec. IIIII In Sec. IIVI we present the phase 
diagram and discuss its implications. Sec. |V]is devoted 
to the multifractal analysis of the atomic wavefunction, 
from which the proposed phase diagram is verified. The 
diffusion dynamics is studied in Sec. I VII for a initially 
localized Gaussian wavepacket. Finally we conclude in 

Sec. ivni 



II. MODEL AND HAMILTONIAN 

We consider the following model depicted in Fig. [1] 
cold atomic gas with internal spin states |t) and jj,) 
confined in a spin-independent ID quasiperiodic optical 
lattice potential siEfn sin^ {kix) + S2-E_ri sin^ {k2X + 4>) 
along the a;-direction, which is formed by combining 
two incommensurate optical lattice together Here 
ki = 2n/Xi is the lattice wavenumber, Si is the height of 
the lattice in unit of the recoil energy Em — h? /2m\\ 
and <j) is the relative phase between the two standing- wave 
modes (without loss of generality, we assume = in the 
following discussion) . It is assumed that the depth of the 
lattice with wavevector ki is deep enough to serve as the 
tight-binding primary lattice. In the meanwhile, a pair 
of counter propagating Raman beams couple the atomic 
states It, kx = q) and Ij, kx = q + 2kji), which creats the 
effective SO coupling |13j . 

In the basis composed of atomic pseudo-spin states 
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FIG. 1: (Color online) Schematic diagram showing the system 
under consideration. 



{It) : the single-particle Hamiltonian reads 



2m 



SiEjii sin^ {kix) + S2ER1 siv? {k2x) 



n 





-2ikn_x 







with Q, the effective Raman coupling strength. Here 
we have assumed that the Raman two-photon detun- 
ing is 0. Then introduce the dressed pseudospin states 
{It)' = It) e-^'=«^, It)' = |t)e*'=«^}, which is equivalent 
to performing a pseudo-spin rotation with the operator 
R{kii) = exp (— ifcflXiTz), and perform a global pseudo- 



spin rotation (t^ 



ay, ay — > ax, ax — >■ a^ I13|, |1 



In 



the new basis, Hamiltonian h can then be rewritten as 

h = hso + V 



n 



2m 



^(T2 + si£^_Ri sin^ (kix) 



S2ER1 sin^ {k2x) , 



(1) 



in which hso describes single-particle motion in 
the presence of the effective SO coupling, which 



is embodied in the vector potential A 



-mXa,„ 



(A — hkji/m characterize SOC strength) and the effec- 
tive Zeeman field r2/2. 

h effectively describes a spin- 1/2 particle moving in a 
ID quasiperiodic potential and subject to both the Zee- 
man field and equal Rashba-Dresselhaus SO coupling, 
which can be used to simulate the corresponding con- 
densed matter system such as the motion of electron in 
one dimentional semiconductor nanowire with disorder 
and SO interactions. 



III. TIGHT-BINDING APPROXIMATION 

In the language of quantum field theory, the total 
Hamiltonian describing our system reads 



Hix)= / dx*^ (x) h (x) if (x) , 



(2) 



with 



the atomic field operators. 



In the tight-binding limit, the atomic field opera- 
tor can be expanded as "i! (x) — '^jWj{x)cj, where 
■"^i i^) = w [x — Xj) is the Wannier state of the primary 
lattice at the j-th site and Cj — {cj^,Cji) are annihila- 
tion operators. By considering that the tunneling takes 
place between nearest neighbour sites and retaining only 
the onsite contribution of the secondary lattice, one can 
achieve the following description of ([2]) (with the energy 
measured in units of Em and length scaled in units of 



H = J2 - (c]rcj+i + H.c.^ - A cos (27r/3j) c]cj + ^c]a,Cj 



< — J c] (cos TT-f — idy sin7r7) Cj+i -I- H.c. — A cos (27r/3j) c^-Cj + —CjazCj 



where the tunneling amplitude T — Jexp is obtained through Peierls subsititution [l^, [23| , which was 



also used in recent works to study the impact of SO coupling on two-dimensional Bose-Hubbard model |21|-|2 
J Adl = A [xj — Xjjf^i) = tvK^ay is the integral of the vector potential along the hopping path with 7 = kn/ki. J is 
the tunneling amplitude without SO coupling, /3 = k2/ki. J and A can be calculated as 



J — — J dxwj+i (x) 

S2P^ 



dx^ 



si sm X 



Wj (x) 



A = 



dxcos {2l3x) \w (x)l^ 
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By writing 



Si a V''j,o-c| the stationary Schrodinger equation H\i\)) = E \tp) lead to 



-J cos(7r7) (V'i+i.t + V'j-i,t) - sin(7r7) (i/'j+i,; - V'i-i,;) - A cos (27r/3j) V'j,t + ^V'j.t = ^V'j.t > (3a) 
-J cos (7r7) (V^j+i,^ + -0^-1,;) + ./sin (7r7) - ipj-i.t) ~ A cos (27r^j) -0^4 - ^0j,_^ = i;^'^,^ . (3b) 



r 



The second terms on the left hand side of Eqs. (O, 
which are proportional to Jsin(7r7), represent the spin- 
flipping tunneling which arises from the effective SO in- 
teraction. In the absence of SO coupling (7 = 0, 17 = 0), 
spin-^ and 4- components are decoupled and we have 

- J (V'j+i + 0j_i) - Acos (27r/3j) 0^ = Et^j, (4) 

which represents the typical Harper equation [24|] or the 
Aubry- Andre model [23] ■ Eq. ^ satisfies Aubry duality, 
as can be seen by performing the transformation ipj = 



- y (^m+l + "^m-l) - 2 J COS (27r^m) 0„ = i^l/'m, (5) 

Eqs. dU and (O are identical at A/ J = 2. Since the 



transformation made above represents the typical Fourier 
transform which transforms localized states to extended 
states and vice versa, then the critical point A/J = 2 
is identified as the transition point between the local- 
ized states and extended states, i.e., all the single-particle 
states are extended when A/J < 2 and localized when 
A/J > 2. 

The properties of the Aubry-Andre model, as repre- 
sented by Eq. ^ , have been theoretically studied [1^ [l^l 
and it can be implemented in systems of Bloch electrons 
(20} and cold atoms [5] . We have also studied this model 
with A dressed by a cavity field through nonlinear feed- 
back 1]. 



A similar transformation ipj,cr = Co- Em 0m, o-^ 
(e^ = 1, ej, = i) made to Eqs. @ will lead to 



y (V'm+i,t + 4'm-i,t) + "^J sin (n-y) sin (27r/3m) -0^4 - 2 J cos {■wy) cos (27r/3m) 0^^^ -|- — 0m,t = /^0m,t i (6a) 
y (^0m+i4 + 0m-i,4,J + 2 J sin (777) sin (27r^m) i/'m,t ~ 2 J cos (7r7) cos (27r/3m) 0™,; - — 0m4 = -E0m4 ■ (6b) 



A comparison between Eqs. (|3]) and Eqs. ([6]) shows that 
the presence of the spin-flipping tunneling terms breaks 
the duality. This distinguishes our current work from 
that reported in Ref. [ll| , in which the authors studied a 
system of two-dimensional (2D) electrons on a square lat- 
tice subject to Rashba spin-orbit coupling and immersed 
in a perpendicular uniform magnetic field. In this system, 
it has been shown [ll| that a gneralized Aubry duality 
is preserved when tunneling along the two perpendicular 
lattice directions are exchanged. Such an operation is 
not available in our system as ours is an intrinsically ID 
model. 



IV. PHASE DIAGRAM ANALYSIS 



Here we follow the method in [Tl[ to map the phase 
diagram using a quantity called the total width of all the 
energy bands B, which had been proved to be useful in 
investigating phase transition in a quasiperiodic system 
[Til [261 1 2 7.] . In order to observe its property, as people 
usually do, we first choose the optical lattice wavelength 
ratio /3 to be /3„ = Fn/Fn+i = p/q, in which F„ is the 
n-th Fibonacci number defined by the recursion relation 
F„+i = Fn + Fn-i with Fo = Fi = 1. When 71 00, 
Pn (\/5 — 1) /2, which is the inverse of the golden 
ratio. 



Due to the lack of duality in the current model, it is 
not clear whether there exists a phase transition between 
the localized and extended states in the present system. 
In addition, what effect will SO coupling take? Will it 
enhance the tendency to localization or delocalization? 
We focus on these problems in the following discussion. 



Since p and q are integers prime to each other, 
the system is periodic with the period q. Under 
the periodic boundary condition, according to Bloch's 
theorem, ipi+q^a = e**'^^0i^CT- Eqs. ([3]) then re- 
duce to an eigenvalue problem Hip — Eip with 

= ('01,t> 014' "^24, 024, • ■ •, 0'g-l,t' 09-14' V'g,!' ^q,-l) 

and the 2q x 2q matrix H takes the following form 
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n = 



(Hi L ■■■ 

it H2 L 

H3 L 

••• 



e-*'=-«Lt \ 




L 







L 



with 



_ , -Acos(2^/3j+0) + f 





-A cos (27r/3j 



and L = 



The Hermite 



— Jcos(7r7) — Jsin(7r7) 
Jsin(7r7) — Jcos(7r7) 
matrix "H can be diagonahzed with 2q real eigenvalues 
Ei (kx), which form 2q energy bands as the function of 
kx in the first Brillouin zone q\kx \ < tt. 

In the absence of SO coupling, the energy bands are 
degenerate for spin-^ and I, with the band edges locate 
at fca; = and k^ = ±7r/q, so B can be calculated hy B = 

Y^=i (0) ~ (i""/?)!- It was first demonstrated in 
[2a | that for extended states with A/J < 2, B approaches 
a finite value for q — )• 00; while for localized states A/J> 
2, B rapidly tends to zero as q ^ 00; at the critical point 
A/ J = 2, B ^ q~^. In this manner, the critical value 
A = Ac signaling AL transition can be determined by 
observing the property of i? as a function of the period 



Q. 



□ 
□ 

□ 

0.4 □ 



0.2 



critical 



o 

o 



0.0 - 



extend 



1.5 



2.0 

A 



° localize 



2.5 



FIG. 2: (Color online) Phase diagram for Anderson localization 
of SO coupled BEC in ID quasiperiodic lattice in the parameter 
space A-f2. 7 = 0.7, A and f2 are estimated in units of J. 

Now taking SO coupling into account, we anticipate 
that the phase diagram is composed of three different 
phases: (i) Extended phase in which the energy spectrum 
is purely continuous and all the eigenstates are extended; 
(ii) Localized phase is characterized by purely dense point 
and all the wavefunctions are localized; (iii) The energy 



spectrum is mixed in the critical phase with extended 
eigenstates coexist with localized ones. Using the method 
described above, Ac is determined as a function of fl, 
which is indicated in the phase diagram of Fig. [5] by 
the line separating the regions signaling localized phase 
and critical phase. Calculation is performed with the 
parameter of 7 — 0.7, which is used throughout the paper 
and can be experimentally realized for *^Rb atoms by 
adjusting the angle between Raman beams [l3|. At = 
0, AL transition occurs at A/ J = 2, reminiscent of the 
situation without SO coupling. This can be understood 
from Eq. ([l} by that SO interaction can be removed from 
the Hamiltonian through a unitary transformation with 
the operator S — exp [—ixay/2S) when = 0. Examples 
of data are shown in Fig. ^&). At A/ J = 2.02, B tends 
to zero for Vl < Vt (Ac) and B tends to a finite value for 
17 > 17 (Ac). 

Due to that duality is broken by the SO coupling here, 
the boundary between extended phase and critical phase 
is not related to that separates localized phase and crit- 
ical phase, which is different from [ll|. The extended 
phase and critical phase cannot be differentiated by ex- 
aming the properties of i3, since the energy spectrum 
contain absolutely continuous parts in both these two 
phases, which leads S to a finite value in the quasiperi- 
odic limit, as shown in Fig. I^Kb). 
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FIG. 3: Total bandwidth versus period q. At g — >■ 00 the system 
becomes quasiperiodic. The parameters are specified in the figure. 

The localization property of the atomic wavefunction 
can be characterized with the inverse participation ratio 
(IPR) which is defined as 

N 



in which N is the number of lattice sites, '0j,t(4.) ^re so- 
lutions of Eqs. ([3]) and fulfil the normalization condition 

Ylij (lV'j,tl^ + = 1- IPR. reflects the inverse of the 

number of the lattice sites being occupied by the atoms. 
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For extended states, — >■ 1/N and approach for 
large N. While for localized states, IPR tends to a finite 
value and a larger value of means that the atoms 
are more localized in space. 




FIG. 4: Eigenenergies obtained from numerical diagonalization 
of Hamiltonian matrix, as a function of A/ J. Calculations are 
performed under periodic condition with /3 — 610/987, A'' = 987, 
7 = 0.7 and (a) = 0.1; (b) = 1. 



V. MULTIFRACTAL ANALYSIS OF 
WAVEFUNCTIONS 

To check the proposed phase diagram, we investi- 
gate the scaling property of the wavefunctions using the 
method of multrifractal analysis described in Take 



the period of the lattice to be Fibonacci number F„, 
from the wavefunctions {ipj.a} we have the probability 

Pj ~ iV'j.tl^ + l'^/'j-J.I^' which is normalized as ^^ZiPj — 1- 
The scaling index a for pj is defined as pj = F~". We 
then assume that the number of sites satisfying the scal- 
ing is propotional to fI"^°'\ / (a) can be calculated as 
/ (a) = lim„^oo /n (a). 

The localization properties of wavefunctions are char- 
acterized by / (a) in the following manner: For extended 
wavefunctions, all the lattice satisfy pj ~ F~'^ , so / {a) 
is fixed at / (a = 1) = 1. On the other hand, a local- 
ized wavefunction has a nonvanishing probability only 
on a finite number of sites. These sites have a = 
[/ (0) — 0] and other sites with probability zero have 
a = oo [/ (oo) ~ 1]. For critical wavefunctions, a has 
a distribution, which means that / [a) is a smooth func- 
tion defined on a finite interval [amin.amax]- 

Numerically we calculate /„ (a) for Fibonacci indices n 
and extrapolate them to n — >■ c». One can then discrimi- 
nate extended, localized and critical wavefunction by ex- 
aming the minimum value of a in the following manner 

extended wavefunction amin = 1, 

critical wavefunction amin 7^ 0, 1, 
localized wavefunction amin = 0. 



We use IPR to determine the boundary separating the 
extended phase and critical phase, which is identified by 
the turning point of IPR as a function of A /J, as those 
had been done in many previous works [^,[l3,[H,|33]. The 
calculation is performed with /3 = F14/F15 = 610/987 
and N = _Fi5 = 987 under periodic boundary condition. 
This give the extend-critical phase boundary shown in 
Fig. [21 which indicate that with the increase of the Rabi 
frequency O the system is more likely to start become lo- 
calized. In order to understand this, we plot the energy 
spectrum as a function of A/ J. When J7 is relatively 
small, the spectrum shown in Fig. Hl^a) possesses similar 
properties as that in the absence of SO coupling: Along 
with the increase of A/ J, two major gaps devide the spec- 
trum into three parts, each of which in turn devides into 
three smaller parts, and so on. This is because the value 
of 1//3 lies between 2 and 3. The spectrum of localized 
states is then characterized by the presence of an infinite 
number of gaps and bands. The effective Zeeman term 
which is propotional to 51, in combination with SO cou- 
pling and the lattice structure, take the effect of opening 
gaps between different energy bands, as shown in Fig. 
IDJb) . So the critical value Ac takes a smaller value with 
the increase of VI. 
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FIG. 5: amin versus 1/n for the wavefunctions of the lowest band 
O (« = 1) and the centre band U [i ^ F„). (a) Q/J = 0.2, 
A/J = 1.5; (b) n/J = 0.4, A/J = 1.5; (c) Q/J = 2.5, A/J = 
0.1 corresponds to extend phase, critical phase and localize phase, 
respectively. 

Q^min is calculated for example wavefunctions and the 
results are shown in Fig. [5] The wavefunction of the 
lowest band is denoted by i = 1, while that at the centre 
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of the energy spectrum is denoted by i = F„. First, for 
fi/J = 0.2, A/ J = 1.5 at which the system is in the ex- 
tended phase according to the phase diagram Fig. [21 ckmin 
extrapolates to 1 for both i — 1 and i — Fn, as shown in 
Fig. [Sja), indicating that the energy spectrum is purely 
continuous and all the wavefunctions are extended. The 
point ri/J = 0.4, A/ J = 1.5 corresponds to the critical 
phase in Fig. [51 and in Fig. [5]Jb) one can find out that 
ttmin extrapolates to for i = 1 and amin extrapolates 
to 1 for i = J^„. This suggests that the wavefunction at 
the edge of the energy spectrum is localized while that 
at the centre is extended, which indicates the existance 
of mobility edges. 

The appearance of mobility edges here can be under- 
stood as the result of the breaking of original self-duality 
via SO interaction, which can also be aroused by other 
effects such as hopping beyond neighbouring lattice sites 
[l^lso'l. We would like to note that, even if the duality is 
preserved, SO coupling can also lead to the appearance 
of critical phase and mobility edges, as those had been 
demonstrated in [TTj . 

VI. DIFFUSION DYNAMICS 



is intimately related to the localization properties of the 
system: 

(i) in the extended phase of A/ J < 2, the wavcpacket 
will experience ballistic expansion with 77 = 1; 

(ii) at the critical point of A/ J = 2 the wavepacket 
subject to subdiffusion with rj ^ 0.5; 

(iii) the wavepacket is localized when A/J > 2 and 
ry = 0. 



without SO interaction with SO interaction 




log,o(Jt) log^„(Jt) 



In realistic experiment, localization properties can be 
investigated by loading the SO-coupled BEC into the 
quasiperiodic potential and observing its transportation 
along the lattice (H . The equations-of-motion associated 
with Eqns. ^ are 



= -Je- 



-Acos(27r;3j) + 



(7) 



where = (V'j.t' V'i.J.) ■ We study the diffusion of ul- 
tracold atomic gas in quasiperiodic optical lattice with 
Eq. ([7]) by taking the initial atomic wavefunction to be 
a localized Gaussian wavepacket with width a 



(< = 0) = (2a^/7) 



in which we take a = 5 in the following calculation. Here 
we assume that the atomic wave packets initially lie in the 
centre at j = with the system size of 2000 lattice sites. 
The numerical simulation is performed with vanishing 
boundary condition and during the time evolution the 
atomic wavepacket never reaches the boundaries so that 
the effect of boundary condition does not appear. 

To measure the localization, we use the width of the 
wave-packet defined as 



w 



(Ax) 




In the absence of SO coupling, the time evolution of w (t) 
can be parametrized as w {t) ^ @,|3lii a-nd its property 



FIG. 6: Time evolution of the wavepacket width w (t). (a) without 
SO interaction; (b) in the presence of SO interaction with 7 = 0.7 
and n/J = 0.2. 

The above time-evolution properties are demonstrated 
in Fig. [H^a). In Fig. [S^b) we present the results in 
the presence of SO interaction. One can find out that, 
for A/J = 2.1, which corresponds to the critical phase 
shown in Fig. [21 the wavepacket still subdiffusc with time 
evolution. In addition, the time evolution of wavepacket 
at sample time are shown in Fig. [7l At A/J =1.5 
correspond to the system in the extended phase, the 
wavepacket rapidly diffuses and almost all the lattice 
sites become populated. While at A/J — 2.5 for the 
localized phase of the system, there are no diffusion be- 
cause in this regime the initial Gaussian wavepacket can 
be decomposed into the superposition of several single- 
particle localized eigenstates. For the system in the crit- 
ical phase at A/J = 2.1, the wavepacket diffusion is ac- 
companied with solitonlike structures in the centre and 
spreading sideband, which reflects that extended and lo- 
calized eigenstates coexist in the system. 

Besides that, the nature of localized states can also 
be extracted from the momentum distribution of the 
atom stationary states. This is because a more localized 
atomic wavefunction corresponds to a wider momentum 
distribution via Fourier transformation. The momentum 
distribution can be measured by turning off the Raman 
lasers, releasing the atoms from the lattice and perform- 
ing time-of-flight imaging. Since our above discussions 
are for the dressed spin states {|t, k)^ , \l, k)^}, their mo- 
mentum distribution can be mapped from that of the 
bare spins |t, k + kn) and \\., k — kn) in the following 
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Jt=50 Jt=500 Jt=1000 




lattice sites lattice sites lattice sites 



FIG. 7: Time evolution of wavepackets |i/'j,tl correspond to Fig. \Elh) at specific times. 



manner 

It: = ^ It, k + kR) + i;, k-kn), 

It' fc>d = ^ It, k + k„) + k - kli) . 

VII. CONCLUSION 



alization of the system using cold atomic gas trapped 
in a quasiperiodic optical lattice potential and external 
laser fields. Since the two ingredients of our proposed 
scheme, the quasiperiodic optical lattice potential [5| and 
SO coupling [l^, [ll, [l^, had already been achieved for 
cold atoms, it is expected that the localization properties 
discussed in this work can be readily observed in experi- 
ment. 



In conclusion, we have studied the system of a SO- 
coupled spin-1/2 particle moving in a one-dimensional 
quasiperiodic potential. We mapped out the system 
phase diagram in the tight-binding regime and accord- 
ingly discussed the localization properties. In the ab- 
sence of SO interaction the system can be mapped into 
the AA model and self-dual if A/ J = 2, SO interac- 
tion breaks the duality and leads to the appearance of 
critical phase, in which the extened and localized states 
coexist in the energy spectra. We also verified the phase 
diagram via multifractal analysis of the wavefunctions 
and diffusion dynamics of a initially localized Gaussian 
wavepacket. Experimental detection of localization prop- 
erties are discussed. We proposed an experimental re- 
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